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A weak-localization effect has been studied in macroscopically inhomogeneous 2D system. It is 
shown, that although the real phase breaking length tends to infinity when the temperature tends 
to zero, such a system can reveal a saturated behavior of the temperature dependence of that 
parameter, which is obtained from the standard analysis of the negative magnetoresistance and 
usually identified by experimentalists with the phase braking length. 
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The numerous experimental as well as theoretical pa- 
pers are devoted to the weak localization in semiconduc- 
tor structures and metal films. The problem was most 
intensively studied in the eightieth years. The theory 
of weak localization, which adequately describes exper- 
imental data, was developed. It gave simple analyti- 

' cal expressions for quantum correction to conductivity 
which allowed to determine a phase breaking time in real 
electron systems experimentally. The comprehensive re- 
view of the status of the problem at that time is given in 

■ Ref. 1,2. 

In recent years the interest to this problem reappeared. 
One of the reason is that the new experimental results 
have been obtained. One of them is the saturation of 
the temperature dependence of the phase breaking time 
T^p.^ These papers cause a storm discussion in the liter- 
ature (see, e.g., Refs. 5,6,7 and references therein) and 
\ stimulate a new flux of the papers concerning the weak 
localization. 

A standard fitting procedure is used to analyze the ex- 
perimental data and determine the phase breaking time 
practically in all cases: experimental magnetic-field de- 
pendencies of magnetoresistance are fitted to the the- 
\ oretical ones,^'^ the phase breaking time is the fitting 
parameter. 

Another approach to examination of the negative mag- 
netoresistance due to weak localization has been pre- 
\ sented in Ref. 10. It is based on a quasi-classical treat- 
ment of the problem^^'^^'^'^ and an analysis of the statis- 
tics of closed paths of a classic particle moving with 
scattering over 2D plane. This method has been ap- 
plied to study the weak localization in InGaAs / GaAs het- 
erostructures with singie^^ and double^^ quantum wells. 

The foregoing relates to homogeneous systems reveal- 
ing only weak disorder on microscopic length scale. Of 
special interest are inhomogeneous systems, e.g., gran- 
ular CuO films^^ and percolating gold films. ^^'^^'^^ For 
interpretation of the experimental results the concept of 
anomalous diffusion is widely used. It should be men- 
tioned that the transition from the weak localization 
to strong localization regime, when the system param- 
eters or external conditions are varied, can also cross the 
macroscopically inhomogeneous state in originally homo- 



geneous 2D systems. 

In the present paper we analyze the weak-localization 
negative magnetoresistance in macroscopically inhomo- 
geneous systems using approach developed in Ref. 10. 
The purpose is to show how the parameters, determined 
from the magnetoresistance experiments, match their 
real values in the systems with different inhomogeneity. 
Moreover, we show that this approach allows experimen- 
tally to distinguish what geometry of inhomogeneity is 
dominant in the system under investigation. As in Ref. 
10 the results presented have been obtained from the nu- 
merical simulation of classical motion of particle. 

We suppose that a macroscopically inhomogeneous 2D 
system consists of a number of puddles which are con- 
nected one with other by means of channels (Fig. 1(a)). 
The transport through the puddles and channels is dif- 
fusive, i.e. their dimensions are much greater than the 
mean free path of electrons. In this case the quasi- 
classical approach to consideration of the problem can 
be applied. We focus our attention on the case when the 
system is in inhomogeneous regime, i.e, when L,^ < ^p,^" 
where S,p is the percolation correlation length, = 
^JDt\p, the diffusion constant is given by £> = v\tI2 
with vf as the Fermi velocity and r representing the 
transport elastic mean- free time. Furthermore, we have 
idealized the situation and supposed that all the chan- 
nels and puddles are identical, and our model system is a 
series-parallel connection of the elements, each of them is 
a channel and two half-puddles connected (Fig. l(b,c)). 
We reduce thus the problem on the infinite inhomoge- 
neous 2D system to the problem on the weak localization 
in one element, which is a diffusive constriction connect- 
ing two diffusive puddles. Because the channel and two 
paddles are connected in series, the correction to the con- 
ductance of the element is 

^a^^—f^^^ + ^^^V (1) 
Rc + Rp \i?c + Rp Uc Rc + Rp cTp ) 

where i?c, Rp are the resistances, CTc, fp and (5ctc, iop are 
the conductivity and correction to the conductivity, re- 
spectively. Indexes c andp, hereafter, refer to the channel 
and puddle, respectively. 

The quantities 8ac and are determined by 

the statistics of closed paths and can be found as 
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FIG. 1: (a) Sketch of inhomogeneous 2D system. The key 
parts determining the conduction are enclosed in rectangles. 

Non-conducting areas arc white. Points represent scatterers. 
(b), (c) The models of long and square channels, respectively, 
used in simulation. Polygonal lines show particle trajectories. 
Shadowing in (c) shows the area with smaller value of the 
Fermi quasimomentum. 

5c7c,p = -2TTllpGoWc,p (2) 

where I denotes the mean-free path within the channel 
(Ip) or puddle (Ip), Go = e"^ / {2^"^ h) , W stands for the 
classical quasi-probability density for an electron to re- 
turn to the start point within the channel (Wc) or puddle 
(Wp). With (7 = irkplGo, where kp is the Fermi quasi- 
momentum, we find 

where K,, and Kp are the length to width ratio of channel 
and puddle, respectively. Using formalism of Ref. 10, we 
can write for qusi-probabilities Wc, Wp in the presence of 
external magnetic filed and inelastic scattering processes, 
destroying the phase coherence: 

>V.(«, = £.S«',5,exp(-|)».(H^), 

(4) 

where Wi{S) are the area distribution functions of closed 

paths, starting within the channel {i = c) or puddle 
{i = p), defined in such a way that Wi{S)dS gives 
the classical probability density of return to the starting 
point with area enclosed S from the interval {S, S + dS), 
lip = VpT^, Li is the average length of closed paths with 
a given area, calculated by appropriate manner [Eq. (6) 
in Ref. 10], $o is the elementary flux quantum. Expres- 
sions (3) and (4) allow to obtain the magnetic field and 
temperature (through Z^) dependence of conductance of 
our inhomogeneous model system. 

Thus, the relative contributions of channel and puddles 
to the interference correction are determined not only by 
the statistics of closed paths, but geometry of channel 



and puddles and the values of the Fermi quasimomentum 
in them. 

In the present paper we consider the situation when the 

channel, which connects the square puddles (Kp = 1), 
mainly determines the conductance of whole 2D system 
and interference correction to it in actual range of pa- 
rameters, i.e., the first term in Eq. (3) is greater than 
the second one. Two opposite cases are analyzed: (i) 
fc^ = fc^, /Tc » 1 (long channel); (ii) fc^ < fc^ , = 1 
(square channel). In the second case the lower value 
of the quasimomentum within the channel leads to aris- 
ing of quasimomentum dependent reflection of electrons 
from the opened ends of channel. Those electrons, which 
attempt to enter the channel at an acute angle to the 
channel-puddle border (9 < 9'), arc reflected and return 
to the puddle. All the electrons escaping the channel 
do this freely without any reflection. The value of 6' is 
determined by hp to kp ratio. The reflection can be spec- 
ular, if the border is smooth, or difi^usive, if the border 
is ragged. For simplicity we suppose that the border is 
straight line and reflection is specular. 

Let us describe the simulation procedure. The element 
is represented as a lattice which contour corresponds to 
the geometry of element. The scatterers with given cross- 
section are placed in a part of lattice sites with the use 
of a random number generator. A particle is launched 
from some random point, then it moves with a constant 
velocity along straight lines, which happen to be termi- 
nated by collisions with the scatterers. After collision it 
changes the motion direction. If the particle collides with 
the walls of element, it is specularly reflected. If the par- 
ticle passes near the starting point at the distance less 
than d/2 (where is a prescribed value, which is small 
enough), the path is perceived as being closed. Its length 
and enclosed algebraic area are calculated and kept in 
memory. The particle walks over the element until it 
reaches one of the channel belonging to another element. 
As this happens one believes that the particle has left to 
infinity and will not return. A new start point is chosen 
and all is repeated. 

The parameters used in simulation are the following: 
d = 5 [hereafter we measure length and area in units of 
lattice parameter and its square, respectively]; the num- 
ber of launches Ig = 10^.. 10^; the length and width of 
channel are 4000 and 400 for long channel, 400 and 400 
for square channel. The dimension of each puddle is 
6000 X 6000 for all cases. The scattering is isotropic, 
the cross-section of the scatterers is equal to 7. The den- 
sity of scatterers is such that the mean-free path is about 
40 for both cases. For the case of square channel we use 
{kl^/k'pf = 10, that gives the value of 9' about 70 de- 
grees. For illustration, let us set the lattice parameter 
equal to 5 A. Then our model provides an example of 2D 
system with local concentration of scatterers 1.5 x 10^^ 
cm-2, ; = 200 A and Btr ^ 0.8 T, where Btr = hc/{2eP) 
is so called transport magnetic field. 

We first turn to the area distribution function W{S) 
(Fig. 2). The area range we are interested in is 5 > Z^, 
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FIG. 2: (a) The area distribution functions Wi{S) (a) and 
3>(S') function (b). Symbols are the simulation results for 
the case of long (o), square (□) channels, and for the paths 
starting within the puddle (a). Only the positive area range 
is presented. Results obtained for negative algebraic areas are 
identical. 



when the approach to analysis of the negative magnetore- 
sistance suggested in Rcfs. 10 is applied. As is seen the 
<I'-versus-S' curve shows more rapid than decreasing 
in wide area range for both cases considered. 

We turn now to consideration of the interference quan- 
tum correction to the conductivity. It has been found 
using Eq. (3) and discrete form of Eq. (4) [see Ref. 10 for 
details] : 



1 



Ldl 



El bSk 



exp -7- 



c,p. (6) 



Here, summation runs over all closed paths, Sk, Ik are 
the area enclosed and the length of fcth closed path, re- 
spectively, b is magnetic field measured in units of Btr. 
Deriving Eq. (6) we have supposed that r<^ is the same 
within the channel and puddles. In this case Wi are in- 
dependent of the Fermi velocity vp- Figure 3(a) shows 
the simulation results for 6a {b = 0), obtained for different 
I /lip values, and, for comparison, the results of theoretical 
calculation obtained through the well-known formulas^'^ 



because just the long closed paths determine the low 
field negative magnetoresistance at low temperatures. 
As is seen, 4:TTl^Wp{S) mostly follows dependence. 
This is in agreement with the result of diffusion the- 
ory for infinite 2D system, '^'^'^'^ which gives ATrl'^W{S) = 
tanh(7r5'/^2) ~ S'^ for S > 

For the channels, the dependence A'kI'^W{S) looks 
more complicated. For the long channels it is close to 
S^^-^-law within the area range 1 x lO"' — 5 x lO''. It is 
precisely this behavior that is theoretically predicted for 
diffusive motion over infinitely long, narrow strip. A par- 
ticular interest is drastic decrease of the area distribution 
function evident for S > 10^. The origin is that the par- 
ticle, moving along a long trajectory, reaches one of the 
ends of long channel, escapes the channel, and carries on 
the motion mainly within the puddles. 

As for square channel, the 47r/^Wc-versus-5' plot shows 
the behavior close to up to 5 ~ 10*, and reveals 
strong decrease at large areas as in previous case. This 
decrease is caused by the fact that a particle crosses the 
channel/puddle border freely only in one direction: when 
it moves from the channel to puddle. It cannot easily 
return back due to reflection from the opened ends of 
channel resulting from the difference of the Fermi quasi- 
momentum in channel and puddle discussed above. 

In contrast to the area distribution functions, the de- 
pendencies L{S) are similar for all three cases. 

In Figure 2 the function 



(5) 



which determines the magnetic field dependence of in- 
terference quantum correction of the element, is plotted. 

Namely this function rather than the area distribution 
function as in homogeneous 2D case will be obtained 



6a 



-Go In 1 



for 2D case, and 



6a = -6.094 GoL-'^^/^ 



(7) 



(8) 



for quasi one-dimensional case with Lc as channel length. 
As is seen the simulation results reveal the behavior close 
to theoretical ones only for considerably large I to l^p ra- 
tios: l/lip > 0.01. Experimentally, this corresponds to 
the high temperatures. With lowering l/l,p, i.e., with 
temperature decrease, the simulation data tend to satu- 
rate as opposed to the theories. It is clear that such a 
behavior is a result of steep fall in ^{S) dependencies at 
large areas (see Fig. 3(b)). 

Let us analyze the magnetic field dependencies of inter- 
ference quantum correction. As an example, in Fig. 3(b) 
we present the dependencies Aa{b) ~ 6a{b) — Sa{0) cal- 
culated with l/l,p = 3.5 X 10~^. We have considered 
these data as experimental ones and fitted them to the 
well-known expression^'^ 



Aa{b) = aGo 



(9) 



where V'(a^) is & digamma function. The parameters a 
and 7 = l/lif, have been used as fitting ones. Just this 
procedure is usually used to determine the value of l,p 
in real 2D samples. Results of the fit are represented 
in Fig. 3(b) by curves. As is seen, Eq. (9) satisfactorily 
describes the simulation data. 

How the value of 7, foimd from the fitting, matches the 
value oil/lip = t/t^p., put in Eq. (6), is shown in Fig. 4(a). 
The fitting procedure gives the value of 7 which is less 
than the value oil/l,p. In the case of square channel. l/l,p- 
dependence of 7 can be described by the power function: 
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FIG. 3: (a) The interference quantum correction to the con- 
ductivity in zero magnetic field as a function of l/l^ obtained 
from the simulation (symbols), calculated from Eq. (7) (up- 
per curve) and Eq. (8) (lower curve), (b) The magnetic field 
dependence of interference quantum correction to the conduc- 
tivity obtained from the simulation with l/l,p = 3.5 x 10~^ 
(symbols). Curves are the best fit to Eq. (9), carried out for 
h < 0.1. Symbolic designations are the same as in Fig. 2. 



7 oc {l/lip)'^^^. As for the case of long channel, the 7- 
veisus-l /lif, plot shows saturation at l/l^ < 10^^. On the 
assumption of ex. T^", a > 0, Figure 4(a) shows qual- 
itative temperature dependencies of the "phase breaking 
time" as it is obtained from the standard fitting proce- 
dure. Thus, the presence of macroscopical inhomogeneity 
in samples can lead to inconsistency between tempera- 
ture dependence of the "phase breaking time" , found ex- 
perimentally, and the real temperature dependence of r^. 
This inconsistency may be qualitative for inhomogeneous 
samples in which the conduction through the long chan- 
nels dominates: the value obtained from the fitting pro- 
cedure can be saturated, whereas the real phase breaking 
time decreases with decrease of temperature and diverges 
at T — !■ 0. Using quota we underline that the "phase 
breaking time" is nothing more than the fitting parame- 
ter here. It should not be always identified with the real 
value of phase breaking time. 

Finally, we present the ^/l^-dependcncc of prefactor a 
(Fig. 4(b)). Absolute value of a has no meaning, because 
its value should depend on size distribution of channels in 
real sample. However, the dependence is meaningful. In 
both cases considered the value of prefactor a strongly 
depends on I to ratio, and the direction of change 
of a with growing l/l^ is determined by the kind of in- 
homogeneity: for long channel, the value of a increases 
with decrease oil/l^, whereas it decreases for the case of 
square channel. 

Note, that the results presented cannot be extrapo- 
lated down to l/ly, = [or T = 0]. The point is that 
iov l/l^ < 2 X 10~^ the value of = ^Jl^ becomes 
greater than the size of the element, which determines 
the correlation length ^p. In this case our system should 




FIG. 4; The fitting parameters 7 (a) and a (b) as function of 
//?^ ratio. Upper line in (a) is 7 = Ijlip-, lower is 7 oc i^ll^pf'^- 
Symbolic designations axe the same as in Fig. 2. 



behave itself as homogeneous one. Thus, it seems rea- 
sonable to say that the ///.^-dependence of 7 must tend 
to 7 oc /~^, and the value of prefactor must be equal to 
unity for L^'^ i^^i^ 10''. 

In conclusion, we have numerically studied the statis- 
tics of closed paths and the negative magnetoresistance 
in macroscopically inhomogeneous 2D systems. We have 
considered the situation when the main contribution to 
the quantum interference correction to the conductiv- 
ity comes from constrictions. Two types of constrictions 
have been considered: long and square channels. We have 
shown that the magnetic field dependence of the negative 
magnetoresistance is well fitted to the Hikami expression 
for both systems. However, the value of the fitting pa- 
rameter 7 does not coincide with the value Ijl^ fed into 
simulation, and, moreover, the 7-versus-Z/Z,p plot is not 
linear. This means that the temperature dependence of 
the "phase breaking time", as it is obtained experimen- 
tally from the standard analysis of the negative mag- 
netoresistance for inhomogeneous 2D system, can differ 
drastically from the real temperature dependence of T^p. 
An indicator, whether the system is inhomogeneous, is 
the strong temperature dependence of the prefactor. The 
Fourier transformation of the magnetic field dependen- 
cies of the negative magnetoresistance taken at differ- 
ent temperatures^"'^'''^^ can be used in order to obtain 
the function <&(S'), which carries direct information about 
the area distribution function of closed paths, and, con- 
sequently, about the character of inhomogeneity of the 
samples, investigated. 

The results presented can be directly used to interpret 
the weak-localization experiments carried out on the ar- 
tificial diffusive cavities and constrictions. To our knowl- 
edge there is deficit of the paper devoted to such objects. 
The main attention was attracted to the ballistic cavi- 
ties and systems with antidot arrays [see, for example, 
Ref. 22]. 
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